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We exploit a recently constructed mapping between quantum circuits and graphs in order to prove 
that circuits corresponding to certain planar graphs can be efficiently simulated classically. The proof 
uses an expression for the Ising model partition function in terms of quadratically signed weight 
enumerators (QWGTs), which are polynomials that arise naturally in an expansion of quantum 
circuits in terms of rotations involving Pauli matrices. We combine this expression with a known 
efficient classical algorithm for the Ising partition function of any planar graph in the absence of 
an external magnetic field, and the Robertson-Seymour theorem from graph theory. We give as an 
example a set of quantum circuits with a small number of non-nearest neighbor gates which admit 
an efficient classical simulation. 



I. INTRODUCTION 

From its early days quantum computing was perceived 
as a means to efficiently simulate physics problem [Tl[2], 
and a host of results have been derived along these lines 
for quantum j3lll[5l[6l[Zl|8liail0l[ni[l2lll3llIl[l5], 
and classical systems [16l[17l[18l[19l|20l|2ll[22ll23l[2ll 
E51 E7] . A natural problem relating quantum com- 
putation and statistical mechanics is to understand for 
which instances quantum computers provide a speedup 
over their classical counterparts for the evaluation of par- 
tition functions [THl [H] • For the Potts model, results ob- 
tained in |23j provide insight into this problem when the 
evaluation is an additive approximation. We provided a 
class of examples for which there is a quantum speedup 
when one seeks an exact evaluation of the Potts partition 
function 

In this work we address the connection between quan- 
tum computing and classical statistical mechanics from 
the opposite perspective. Namely, we seek to find restric- 
tions on the power of quantum computing, by employing 
known results about efficiently simulatable problems in 
statistical mechanics. Specifically, we restrict our atten- 
tion to the Ising model partition function Z , and use a 
mapping between graph instances of the Ising model and 
quantum circuits introduced in |28j , to identify a certain 
class of quantum circuits which have an efficient classical 
simulation. 

Restricted classes of quantum circuits which can be 
efficiently simulated classically have been known since 
the Gottesman-Knill theorem [53]. This theorem states 
that a quantum circuit using only the following elements 
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can be simulated efficiently on a classical computer: 

(1) preparation of qubits in computational basis states, 

(2) quantum gates from the Clifford group (Hadamard, 
controUed-NOT gates, and Pauli gates), and (3) mea- 
surements in the computational basis. Such "stabilizer 
circuits" on n qubits can be be simulated in 0(n log n) 
time using the graph state formalism [3D]. Other early 
results include Ref . [JT] , where the notion of matchgates 
was introduced and the problem of efficiently simulat- 
ing a certain class of quantum circuits was reduced to 
the problem of evaluating the Pfaffian. This was sub- 
sequently shown to correspond to a physical model of 
noninteracting fermions in one dimension, and extended 
to noninteracting fermions with arbitrary pairwise inter- 
actions [32, ji3j 34J (see further generalizations in Refs. 
[SSIISB]), and Lie-algebraic generalized mean-field Hamil- 
tonians [37]. Criteria for efficient classical simulation of 
quantum computation can also be given in terms of up- 
per bounds on the amount of entanglement generated in 
the course of the quantum evolution [35] . 

A result that is more directly related to the one we 
shall present in this work is given in Rcf . [25 , but within 
the measurement-based quantum computation (MQC) 
paradigm. MQC relies on the preparation of a multi- 
qubit entangled resource state known as the cluster state. 
It is known that MQC with access to cluster states is uni- 
versal for quantum computation. Reference |25j consid- 
ers planar code states which are closely related to cluster 
states in that a sequence of Pauli-measurements applied 
to the two-dimensional cluster state can result in a planar 
code state. MQC with planar code states consists of a 
sequence of measurements {Mi, M2, . . . , M„, M} where 
the Mi are one-qubit measurements and M is a final 
measurement done on the remaining qubits in some basis 
which depends on the results of the Mi. Reference [25] 
demonstrates that planar code states are not a sufficient 
resource for universal quantum computation (and can be 
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classically simulated). This fact is attributed to the ex- 
act solvability of the Ising partition function on planar 
graphs. Our results complement the work in [2S], as they 
are provided in terms of the circuit model, and generalize 
to Ising model instances that correspond to graphs which 
are not necessarily subgraphs of a two-dimensional grid. 

Other conceptually related work uses the connection 
between graphs and quantum circuits and the formal- 
ism of tensor network contractions, to show that any 
polynomial-sized quantum circuit of 1- and 2- qubit 
gates, which has log depth and in which the 2-qubit gates 
are restricted to act at bounded range, may be classically 
efficiently simulated [36l [39l . A tensor network is a 
product of tensors associated with vertices of some graph 
G such that every edge of G represents a summation (con- 
traction) over a matching pair of indexes. We also use 
a relationship between quantum circuits and graphs but 
whose construction is quite different [55]. Also, Ref. [H] 
connects matchgates and tensor network contractions to 
notions of efficient simulation. 

Finally, other closely related work was recently re- 
ported in |32] (see also [121I111I1S]), which addresses the 
classical simulatability of quantum circuits. Their results 
use a connection to the partition function of spin mod- 
els, as do we, and they too provide a mapping between 
classical spin models and quantum circuits. Specifically 
pertinent to our work is the fact that they give criteria 
for the simulatability of quantum circuits, using the 2D 
Ising model. That is, circuits consisting of single qubit 
gates of the form 6*^*^=" and nearest-neighbor gates of the 
form e*'^'^^^"'^ are classically efficiently simulable. We 
shall discuss how the nearest-neighbor restrictions can 
be lifted while retaining efficient classical simulatability. 

The structure of this paper is as follows. We begin with 
a brief review of the Ising model in Section [ll] where we 
define the Ising partition function Z. In Section [TTT] we re- 
view quadratically signed weight enumerators (QWGT's) 
and their relationship to quantum circuits, and review 
the relationship between QWGT's and Z. In Section 
IIVI we introduce an ansatz that allows one to associate 
graph instances of the Ising model with circuit instances 
of the quantum circuit model. In this section we derive a 
key result: an explicit connection between the partition 
function for the Ising model on a graph, and a matrix ele- 
ment of the unitary representing a quantum circuit which 
is related to this graph via the graph's incidence matrix 
[Eq. ([35|)]. We then present our main result in Section 
[Vj a theorem on efficiently simulatable quantum circuits. 
The proof depends on the fact that there are algorithms 
for the efficient evaluation of Z for planar instances of 
the Ising model. We also discuss the relation to previous 
work. In Section |VII| we present a discussion and some 
suggestions for future work, including the possibility of 
a quantum algorithm for the additive approximation of 
Z. We conclude in Section |VIII| The Appendix gives 
a review of pertinent concepts from graph theory, and 
additional details, including some proofs. 



II. ISING SPIN MODEL 

We briefly introduce the Ising spin model accompanied 
by some notation and definitions. Let G ~ {E, V) be a 
finite, arbitrary undirected graph with \E\ edges and \V\ 
vertices. In the Ising model each vertex i is occupied 
by a classical spin (Ti = ±1, and each edge £ E 

represents a bond J,ij (interaction energy between spins 
i and j). 

Definition 1 An instance o/ the Ising problem is the 
data A = (G, {Jij}), i.e., A represents a weighted graph. 

The Hamiltonian of the spin system is 

HA{(y) = - J^j^^'^j- (1) 

(''j)e-E 

A spin configuration a = {<7i}[^^ is a particular assign- 
ment of spin values for all \ V\ spins. A bond with Jij > 
is called ferromagnetic, and a bond with Jij < is called 
antiferromagnetic. The probability of the spin configura- 
tion CT in thermal equilibrium for a system in contact with 
a heat reservoir at temperature T, is given by the Gibbs 
distribution: Pa{<j) — -^WA{a), where the Boltzmann 

weight is M^a(ct) = exp[-/3iJA(o-)], /3 = l/kT is the in- 
verse temperature in energy units, k is the Boltzmann 
constant, and Za is the partition function: 

ZA(/3)=^exp[-/3i7A(a)]. (2) 

(Unless there is a risk of confusion we will from now on 
write Z in place of Za(/3) in order to simplify our no- 
tation.) Computation of the partition function is the 
canonical problem of statistical mechanics, since once Z 
is known one can compute all thermodynamic quantities, 
such as the magnetization and heat capacity, by taking 
derivatives of = — /c log Z (the free energy) with respect 
to appropriate thermodynamic variables [48j . 

In this work we restrict our attention to the case 
Jij G { — J, 0, J}, with J > 0, which already gives rise to 
the full complexity of spin glass models and the associ- 
ated computational hardness |46| . For example, with the 
above restriction the problem of computing the partition 
function in the three-dimensional spin-glass is NP-hard 

III. QUADRATICALLY SIGNED WEIGHT 
ENUMERATORS AND THEIR RELATION TO 
THE ISING PARTITION FUNCTION 

Quadratically Signed Weight Enumerators (QWGTs) 
were introduced by Knill and Laflamme in Ref. |49j . 



A problem is called NP-hard if the existence of a polynomial- 
time algorithm for its solution implies the existence of such an 
algorithm for all NP-complete problems. 
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Definition 2 A Quadratically Signed Weight Enumera- 
tor is a bi-variate polynomial of the form 

S{A,B,.,y)^ Y: (3) 

fcekcr A 

where A and B are 0,1-matrices with B of dimension 
n X n and A of dimension m x n. The variable b in the 
summand ranges over 0, 1-column vectors of dimension 
n satisfying Ab = (in the kernel, or nullspace of A), 6* 
is the transpose of b, and \b\ is the Hamming weight of 
b (the number of ones in the vector b). All calculations 
involving A, B or b are done modulo 2. 

Note that the evaluation of a QWGT, given that x 
and y are natural numbers, is in general #P-hard, since 
it includes the evaluation of the weight enumerator poly- 
nomial of a classical linear code [5Dj . 



A. QWGTs from Quantum Circuits 

We shall now review in some detail how QWGT's were 
arrived at in Ref. |49j by considering expansions of quan- 
tum circuits. Let f2 be a quantum circuit formed by 
a temporal ordering of N gates gt, and let [/(fi) = 
nfe=Arfffe — 9n ■ ■ ■ 9i be the corresponding unitary op- 
erator. Note that a universal gate set can be achieved 
by allowing arbitrary rotations about tensor products of 
Pauli operators, i.e., each of the N gates gu can be rep- 
resented as 



where 



cos - /-isin - (76, (4) 



(5) 



1=1 



with n being the number of qubits, such that the Pauli 
matrices are 



1 
1 



1 

1 y ' 

-A / 1 

^11 = = I ^ I , aio = = Q _^ 



rotations with cos(0/2) = 3/5 (i.e., = 2 arcsin(4/5)) 
around operators of weight at most two (the weight of 
CTfe is the number of non-zero pairs of bits in b) (Theorem 
3.3, case (e), of [S2]). Letting 



cos 



a 
1 



a 

7 



so that 7 = \/q;2 + q;'2 ^ we rewrite Eq. M 



as 



5fc = - ("^ - iaVfcJ , 
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(6) 



(7) 



The gate set is still universal if U{Vl) is expressed as 
a product of real gates i.e., if each gate gk contains 
an odd number of cry's, so that iab^ in Eq. ([t]) is a 
real- valued matrix. Following Ref. [IS], we adopt this 
convention, so that from now on bk is a binary vector of 
length 2n, subject to the restriction that the bk can only 
contain an odd number of ll's. Moreover, the gate set is 
still universal if we assume that the orientation (the sign 
of 0) is positive if the number of ay's is lmod(4) and 
negative otherwise \^ . This means that we can replace 
Eq. by 



gk = - {al±ia'abk) : 
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(8) 



with the sign determined by the number of cry 's in atf. 
Then, by defining 



(9) 



where |6|y is the (always odd) number of cry's occurring 
in crf,^ , we may write 



Qk = -{al + a'abj, 
7 



(10) 



which is the desired representation of real- valued gates 

m- 

Now define C to be the block diagonal matrix whose 
blocks consist of 



1 




(11) 



Here b, G {00,01, 10, 11}, b = {bjf^i is a binary vector 
whose length is 2n, and the superscript (i) represents the 
qubit which is operated on by the corresponding Pauli 
matrix. A circuit constructed using gates of the form Q 
may be approximated efficiently to accuracy 0{e/N) with 
polylog(7V/e) overhead using a standard gate set, such as 
controUed-NOT with single-qubit gates, and there is a 
classical algorithm that computes such approximations 
efficiently [321 HI] ■ A universal set of one- and two-qubit 
gates can be obtained from gk's as in Eq. Q from the 



(12) 



Then the property that b has an odd number of ll's is 
given by b^Cb = 1. In addition, we have the multiplica- 
tion rule 



(13) 
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where the addition in the subscript is bit by bit modulo 
2. 

Let H be the (2?i x N) matrix whose columns are the 
bk- 



H = (61 62 



bN). 



(14) 



_ff is a linear size, bijective representation of the quantum 
circuit, where each column represents a gate and every 
pair of rows represents a qubit. 

Definition 3 A matrix H which is constructed according 
to Eq. (14-) is called the "H -matrix representation" of the 



quantum circuit 51. 



Using the rule ( 13 1 we then have the following expan- 
sion US: 



u{n) 



n 9k 

l[ -{al + a'at,) 

k=N ^ 

1_^(_1)Q..«H(«')A^-I.|^^^, (15) 



where 



(16) 



and where the NxN lower-triangular matrix Q is defined 

by 



Q EE Iwtr (iJ'Ci?), 



(17) 



and where a ranges over all binary column vectors of 
length N. (Thus Qaa' is not a matrix element of Q; we 
use this notation merely for convenience.) 

In order to make contact with the partition function 
of the Ising model, we shall be interested in the matrix 
element {0\U{n)\0), where |0) = ^f^jO,), and \0i) is 

the -t-1 eigenvector of ay. For this matrix element to 
be non-zero no qubit can be flipped, i.e., U(fl) cannot 
contain any ax or ay factors. When taking the same 
matrix element of the right-hand side of Eq. ( 15 1 we have 



(O|o'_f/a|0), and similarly, for this to be non-zero ajfa can- 
not have ax or ay factors. This is enforced by summing 
only over those binary vectors a such that CHa — 
Thus: 



{omnm^^ (-i)Q»"ai" 



(18) 



aekor CH 



A glance at the QWGT expression ([s]) reveals a striking 
similarity to the latter matrix element. 



B. Example 

As a simple example meant to illustrate the correspon- 
dence between the if-matrix representation of a quantum 
circuit and the actual operation of the circuit, consider 



1 1 1 

1 

1 1 

1 
1 1 1 
1 1 



H 



This matrix represents a circuit f2 comprising three gates 
(three columns) acting on three qubits (two rows per 
qubit), 

U{n) = 335251, 



with the following unitaries: 

91 = 

92 = 

53 = 



{al-ta'a^z'^a^^^a'^'^) 



{al 
{al 



■ I (1) (2) (3)n 

la a\ 'a^z'ay ) 



■ ' (1) (2) (3)^ 
la CTy (^z "z ) 



The superscripts represent which qubit is being acted 
upon and we have omitted the tensor product symbols. 
The Pauli operators can be read off from the correspond- 
ing column entries in H; thus the entry (1 0)' in the top 
position of the first column of H represents the cr^ Pauli 
matrix in gi, etc. 



C. QWGTs and the Ising Partition Function 

In Ref. ^22] it was shown that the Ising partition func- 
tion can be expressed in terms of a QWGT. Let A be the 
incidence matrix of the graph G = {E, V), i.e.. 



A _ J 1 (w = i and e E) 
^v,i^,3) ]^ else 

Let us associate a binary vector 

W = iwi2,Wi3, ...) 



(19) 



(20) 



of length \E\ with the bond distribution { — ± J}, by 
letting 



1 - J^J/J 



(21) 



so that w specifies whether edge supports a fer- 

romagnetic (wij = 0) or antiferromagnetic {wij — 1) 
bond. Thus we can give an equivalent definition of an 
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instance of the Ising model (recall Definition [TJ as the 
data A = (G, w). 
Let 



A = tanh(/3J), 
and define the \E\ x \E\ matrix 



(22) 



2l^l 


(1 


- A2)l^^l/2 




2l^l 


(1 


- X2)\E\/2 




2l^l 


(1 


- A2)l^l/2 



„ 1 / \ f w on the diagonal 

^ = ^g(^)=|o elsewhere " (^3) 

Writing the instance data as A = (G, w) we then have 
(Theorem 2 of [22 ): 

o\v\ 

^a(a) = E (-1)" 

a^kcr A 

S{A,dg{w),X,l), (24) 
E (-1)"'"AI°I (25) 

atEkcr A 

where a in the sums ranges over all — 1 vectors of length 
\E\ satisfying Aa — 0, where a*Ba — aiWiUi = a ■ 
w (since = or 1) was used in the second equality, 
and where the QWGT definition ^ was used in the last 
equality. 

This establishes the link between QWGTs and the 
Ising model partition function. Because of the similar- 
ity to the matrix element (0|C/(ri)|0), we expect to be 
able to relate the partition function to quantum circuits, 
via QWGTs. We take this up in the next section. 

Definition 4 An even subgraph of a graph G (or equiva- 
lently an Eulerian subgraph) is any subgraph of G whose 
vertices are of even degree. Equivalently, these are paths 
in G which begin and end at the same vertex, and which 
pass through each edge exactly once. 



(26) 



Now, note that the sum in 



5(Adg(w;),A,l)= (-l)""Al''l 



a G kcr A 



is over vectors that are in the kernel (nullspace) of A, 
which here means that only subgraphs having an even 
number of bonds emanating from all vertices are allowed, 
i.e., the sum is taken over all even subgraphs or equiva- 
lently, all Eulerian subgraphs. 

In this work we will sometimes refer to even or Eulerian 
subgraphs as cycles. 



IV. CONNECTING THE ISING MODEL 
PARTITION FUNCTION TO THE QUANTUM 
CIRCUIT MATRIX ELEMENT 



A. A circuit ansatz 

Focusing on the representation of the quantum circuit 
given in Eq. (151 and of the partition function given in 



Eq. ( 25 1 , we begin by asking ourselves if there exists some 



ansatz for the gate set gu such that 



k=N a 



(27) 



where A = tanh(/?J). If such a form were possible, the 
two representations would be closely linked. Indeed, we 
can almost get this form. Let us take as an ansatz 



1 



9k 



(28) 



i.e., the special case of Eq. (10) with A = a' /a, or 

tanh(/3J) = tan(6'/2). (29) 

Note that since the inverse temperature (3 and 
the bond strength J are both positive, Eq. (29 1 
restricts (0/2)mod27r to be in the range (0,7r/2)ij 
(tt, 37r/2), or 0mod47r to be in the range 



R= (0,7r) U (27r,37r) 



(30) 



(the range for which tan(6'/2) > 0). Fortunately, this 
includes the case 9 — 2arcsin(4/5) « 1.85 S R (i.e., 
A — 4/3), which, as noted above, allows a universal set 
of one and two qubit gates to be obtained.^ Thus, we 
have not restricted the generality of the class of quantum 
circuits so far. On the other hand, most 9 Q R do not 
correspond to universal quantum circuits. 



Next, we obtain from Eq. (15 1: 



u{n) 



n 

k=N 



(A2 + 1)^/2 



J2{-l)Q^^X\-\aH^. (31) 



After taking matrix elements (0|-|0), we have, recalling 
Eq. 



1 



(A2 + 1)^/2 



J2 (-1)'3-aI"I. (32) 



aekcr CH 



Comparing Eqs. (25 1 and (32), while using A = a' /a, we 

is to 

A = GH, (33) 



see that a sufficient condition for them to be equal, is to 
identify the incidence matrix A with GH via 



Our goal in this section is to connect the partition func- 
tion Z to {0\U{il)\0). To do so we will use a mapping 
found and described in detail in Ref. [2S] . 



^ These observations were used in Ref. 1281 to show that find- 
ing additive approximations of the signed generating function of 
Eulerian subgraphs over hypergraphs is BQP-complete. 
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where A is a matrix containing twice the number of rows 
of the incidence matrix A, but where each even row is 
a zero row, and each consecutive odd row is equal to a 
row of the original incidence matrix A, and to equate the 
exponents, i.e., find an edge distribution w which solves 



Definition 7 A quantum circuit is "Gw -compatible" 



a ■ w mod 2 = Qaa Va G ker A 
where Q lwtr(iJ*A) [recall Eq. ([itI]. For then 



(34) 



(A2 + 1) 
(1 + A2)^2l^l 



aGkcr A 



(35) 



Equation ( 35 ) is a key result of this paper, as it es- 



tablishes the equivalence between quantum circuits and 
the Ising model, for bond distributions w that satisfy Eq. 
(l34|, and X's that satisfy Eq. (pgl). 



It has two consequences. First, if we are able to deter- 
mine (0|J7(ri)|0), then we are able to determine the par- 
tition function Za(A). Note that estimating (0|J7(ri)|0) 
in general is BQP-complete and thus something one 
could do with a universal quantum computer. Alterna- 
tively, if we had a way of classically computing Za(X), 
then we would be able to classically simulate the quan- 
tum circuit (if it were solving a decision problem) [55] . 
This latter alternative is the one we focus on in this pa- 
per. 



B. Circuit-Ising model compatibility 

The connections we established in the previous sub- 
section between quantum circuits and the partition func- 
tion imply certain restrictions. We flesh these out in the 
present subsection. 

First, since we wish to work only with the physically 
relevant range of positive temperatures and positive J, 
we restrict the gate angles 9 from now on to lie in R. 
Formally: 

Definition 5 A gate angle 9 [Eq. ^] for a gate gk 
[Eq. |jOpy is said to be "\- compatible" if 9 £ R, where 
the range R is defined in Eq. (SC 



Next, we note that Eq. ( 33 1 gives rise to a compatibility 



relation between circuits and graphs: 

Definition 6 A quantum circuit U,, constructed with A- 
compatible angles, is "G -compatible" with a graph G 
if the H -matrix representation of il satisfies Eq. (33), 
where A is the incidence matrix of G, and C is defined 



in Eq. 1711) 



When we take a G-compatible circuit and plug its H- 
matrix into Eq. ( 34 1 we are not guaranteed that there 
exists a solution w. Hence we need an appropriate re- 
striction of the class of G-compatible circuits: 



if it is G-compatible and if the solution set of Eq. (34-) is 
non-empty. 

We need a similar notion for the bond distributions: 

Definition 8 A bond distribution w is "Gfl-compatible" 
with a graph G and circuit f2 if it satisfies Eq. 



34) 



Note that in this last definition the circuit must be 
Gw-compatible, for otherwise we are not guaranteed that 
the solution set of Eq. (34 1 is non-empty. Note further 
that, as these definitions imply, Eqs. (33l-(35) describe 



a connection between quantum circuits and instances of 
the Ising model over given graphs. Namely, any H which 
solves Eq. ( 33 1 is a matrix representation of a circuit f2 



which belongs to a class defined by the incidence matrix 
A of a given graph G. In addition, we can populate the 
vertices of G with weights from the bond distribution w 
provided w is compatible. Thus: 

Definition 9 Let T be any set of graphs for which a so- 
lution to Eq. (34) exists. Then is the set of circuits 
which are Gw- compatible VG G F. 

Definition 10 I{U,rw) is the class of Ising model in- 
stances {A(G,w)}u, whose graph is G Cz T and whose 
bond distributions {w} are Gfl-compatible VG G F. 



Equation ( 34 1 is a system of linear equations over 
GF(2). The number of equations is equal to the number 
of even subgraphs of the given graph or the total num- 
ber of elements in the set ker(G7f), and the number of 
unknowns is equal to the number of edges. However, in 
spite of the fact that the number of elements in ker (Gi/) 
scales exponentially in the number of vertices, it turns 
out that finding a w which solves Eq. (34) can be done 
efficiently [see Eq. ( [38| below] . Let us further stress that 
Eq. (34 1 is only a sufficient condition for the equality of 



Eqs. (25) and (32), and does not capture the whole set of 
possible graph instances that our scheme can handle. We 
define our instances via this condition because it simpli- 
fies the analysis and it allows us to extract information 
about an interesting set of quantum circuits which may 
be classically simulated. We discuss more general suf- 
ficient conditions in Section IVII A| but leave the devel- 



opment of a complete understanding of the actual graph 
instances that our mapping can handle, and in particular 
finding necessary conditions for the equality of Eqs. (25) 
and (32), as a problem for future study. 



V. CIRCUITS CORRESPONDING TO CERTAIN 
PLANAR GRAPHS HAVE AN EFFICIENT 
CLASSICAL SIMULATION 

Let us recap the general idea we have developed so far. 
At the basis of our construction are an inverse tempera- 
ture (3, bond strength J, and a given graph G. We use 
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this graph to first identify a compatible class of quan- 
tum circuits fie (Definition [9| . This class is restricted 
to a subclass of circuits for which there exist so- 

lutions to Eq. (34 1. Such solutions are used to assign 



weights to the graph's edges (a bond distribution), which 
yields a class of Ising model instances compatible with G 
and (Definition 10 1. In other words, we go from the 



unweighted graph to a class of compatible circuits, and 
from there to back to the graph, which is now populated 
by a class of compatible Ising models. Each circuit is 
also parametrized by an angle 9, and when we vary 6 in 
the range R [Eq. (30)] we also vary over /3 and J, via 
tan(6'/2) — tanh(/3jyT^ However, not all values of 9 corre- 
spond to universal circuits. Conversely, not every circuit 
need correspond to a physical (positive) temperature. 

In more detail, we identify the class of quantum cir- 
cuits flc compatible with G (whose incidence matrix is 
A) by solving Eq. ( 33 ) for the matrices H representing 
each (or some) S fie and then find the subset il.Qy 



for which the solution set to Eq. (34 1 is non-empty. We 



then look for a bond distribution w that satisfies Eq. ( 34 ) 



for a given H. Every such w defines an Ising model in- 
stance A{G,w) that is compatible with G and the cor- 
responding n e Qc- We are guaranteed that provided 
such a bond distribution w exists, the partition function 
for the corresponding Ising model is proportional to the 
matrix element (0|[/(il)|0) [Eq. (35 1]. In other words. 



any bond distribution w that satisfies Eq. ( 34 1 induces a 



direct connection between quantum computation and the 
Ising model on a graph with that same bond distribution. 

It is important to emphasize that Eq. ( 34 1 will not 
always have a solution w. Whether or not this is the 
case is entirely determined by the given graph G, since 
G, via its incidence matrix A, determines the class of G- 
compatible circuits ^Iq [i-e., the matrices H that solve 
Eq. (33l], and together they determine a and Qoa that 
go into Eq. (34 1, which w needs to solve. Thus, it makes 



sense to define a class of graphs for which there exists a 
solution w to Eq. ( 34 1 . 



Definition 11 is the set of graphs for which a solution 
to Eq. ( 34 ) exists. 



In order to characterize we require some basic ideas 
from graph theory, such as obstruction sets and down- 
ward closure. These are reviewed in Appendix We 
shall prove: 

Lemma 1 The obstruction set for Q is finite. 

This will come as a consequence of 's downward clo- 
sure and the Robertson-Seymour theorem - Theorem [3] 
This is proved in Appendix [B] We shall also prove there 
that there are no solutions to Eq. ( 34 1 for the graphs K4 



(the complete graph on four vertices) and ^^3,3, where 
K'3.3 is K3 3 with one edge missing (which edge does not 
matter, since upon deletion of another edge one can just 
relabel the edges and nodes and obtain exactly the same 
incidence structure). Formally: 



Lemma 2 The obstruction set for includes 3 and 

The proof is described in Appendix |B] See Fig. |2] in 
Appendix [X] for a pictorial representation of the graphs 
mentioned in Lemma [2j As a consequence we will find 
that all graphs in are planar. 

We shall clarify this conclusion and the previous lemma 
in the next subsection, but given their validity, a Gfl- 
compatible bond distribution w necessarily corresponds 
to a planar graph G. From this it follows that its parti- 
tion function can be efficiently computed classically, and 
hence the corresponding class of quantum circuits, i.e., 
Hqw, also has an efficient classical simulation. Let us 
be precise about what we mean by "classically efficiently 
simulatable" (CES). 

Definition 12 A uniform family Qn — {^i} of n-qubit 
quantum circuits is "classically efficiently simulatable" 
(CES) if the matrix element \ {0\U(fli)\0)\ of each circuit 
in Gn can be obtained to k digits of precision in time 
poly(rt, k) by classical means 155']. 

This definition is a modified version of the one given 
in Ref. [36], which also includes a discussion on how it 
can be weakened. 

When we collect the observations above we arrive at 
an efficient classical test for whether a given quantum 
circuit is CES. This is summarized in Theorem [l] which 
is our main result and the subject of the remainder of 
the paper: 

Theorem 1 (Circuits Corresponding to Certain 
Planar Graphs have an Efficient Classical Sim- 
ulation) 

The class of quantum circuits fieto is CES. Deciding 
whether a given a graph G is in Q can be efficiently de- 
cided. 

The theorem comprises two parts. In the first we char- 
acterize an entire class of CES quantum circuits. The 
proof we offer below is not constructive, i.e., we prove 
that there exists an efficient classical simulation of the 
class of quantum circuits ^qw, and also provide a test 
of non-membership in fie™ for a given quantum cir- 
cuit . In the second part we given an explicit construction 
which decides whether a given graph belongs to the set 
of graphs resulting in CES circuits. To illustrate this 
part, we discuss a class of graphs (which is a subset of 
f2etu) in Section VI for which we can explicitly find the 
GTi-compatible bond distribution. This class is highly re- 
stricted in that the number of even subgraphs only grows 
polynomially in the number of vertices, whereas in gen- 
eral, including restrictions to planar graphs, the number 
of even subgraphs grows exponentially. Nonetheless, the 
class of quantum circuits that one obtains under this re- 
striction is interesting in light of some new results about 
the classical simulatability of quantum circuits [36l |42j . 

For the benefit of the reader we summarize the scheme 
of the first claim of the proof informally. This will also 
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serve to summarize again the mapping between quantum 
circuits and graphs. 

1. Given: any subset F of 0. Every G G T has a 
Gil-compatible bond distribution for some quan- 
tum circuit ri, by assumption. 

2. Take the incidence matrices CH oi the graphs in F 
and transform them into the H-matrix representa- 
tions of the corresponding quantum circuits. The 
following constraint must be respected: Every col- 
umn must have one K-operation and can have at 
most one X-operation. [This constraint comes from 
the fact that CH should be an incidence matrix for 
a graph, where C is the block-diagonal matrix de- 
fined in Eq. ( 12 1 . Without it one has a correspon- 



dence between quantum circuits and hypergraphs 
|28j . Indeed, if the incidence matrix has more than 
two ones per column than one has a hypergraph.] 



3. Thus F corresponds to a set of quantum circuits 
r^rto, i-e., every quantum circuit il e ^Tw is Gw- 
compatible for some G € F. 



4. Show that our mapping from circuits to graphs 
defines a "downward closed set" of graphs, which 
means that we may apply the Robertson-Seymour 
Theorem [S3]. This theorem guarantees that there 
is a finite set of graphs (obstruction set) for which 
we can test whether or not G has any members of 
this set as a graph minor [57j (at most cubic com- 
plexity in the number of quantum gates in il). 

5. Define 0, via this obstruction set, i.e., a graph is a 
member of if it does not have K4 and K^^a = if 3. 3 
with one edge deleted as minors (there may be other 
forbidden minors) . One has a set of circuits which 
correspond to the graphs which have a satisfying 
bond distribution w for equation ( 34 1 . We call the 
corresponding class of quantum circuits ■ (This 
specific obstruction set has been tested with math- 
ematical software. See Appendix [C|) 

6. Due to the fact that these graphs are planar, the 
partition function Z of any graph in Q can be com- 
puted efficiently by a classical computer [50] , 



7. Using equation (35 1, show that knowledge of Z can 



be used to determine the outcome of a quantum 
circuit £ r^rtu for a decision problem. 

8. Conclude: Families of quantum circuits in ileto 
which solve a decision problem can be classically 
simulated. 

Being that F is a subset of Q, any subset of iloto 
is CES. 

Conversely, we have a test for non-membership in the 
set flew- 



1. Input a quantum circuit fl. 

2. Transform f2 into a matrix whose columns represent 
Pauli operations (that are to be exponentiated) and 
every pair of rows are the qubits being acted upon 
as described in Section IIII Al This matrix is called 
H [Eq. ( [I4| ] and is in 1-to-l correspondence with 
n. As above, the following constraint must be re- 
spected: Every column must have one y-operation 
and can have at most one X-operation. 

3. After the above transformation, construct a corre- 
sponding incidence matrix CH of a graph G. 

4. Check the graph G for the minors A''3 3 and K4. 

5. Conclude: If either of these are minors of G then 
reject fi. 



A. Ordering Lemma 

The following lemma allows us to introduce an ordering 
on the elements in Q. 

Lemma S If a graph G is a member of Q, then so is 
G \ Cj or G/cj, i.e., the deletion or contraction of an 
arbitrary edge Cj from a graph in Q is also in &. 

The proof of this lemma is technical and is given in 
Appendix [B] Lemma |3] implies that is a downwardly 
closed set with respect to the minor ordering. Hence we 
can apply the Robertson- Seymour theorem 56J, Theo- 
rem [3] which states that any graph may be tested for 
membership in a given downwardly closed set of graphs 
by just searching the graph for a finite set of minors. The 
complexity of doing this, given knowledge of the minors 
one is looking for, can be shown to be cubic in the num- 
ber of edges. We implemented this to test Eq. (34 1 for 
non- planar solutions, as we describe next. 



B. Equation (34| implies planarity 



Using mathematical software we demonstrated that 
the mapping between graphs and circuits described 
above, with the sufficient condition given by Eq. (34 1, 
cannot be satisfied for and -ftr3,3. That is, if 5 and 
iir3_3 are forbidden minors for 0. The algorithm we im- 
plemented to check this is described in AppendixjC) How- 
ever, a finite graph is planar if and only if it does not have 
Kc, or ^^3,3 as minors (Wagner's theorem; see Appendix 
[A| ). As stated earlier, we thus have 

Lemma 4 All graphs in are planar. 

We remark that we have been able to find examples 
of planar graphs for which there do exist solutions w to 
Eq. (34 1, e.g., ii'2,3- As we explain below, this means that 
includes planar graphs which are not outerplanar. 
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C. Knowledge of the matrix element determines 
output to a decision problem 

The standard way in which a quantum circuit U solves 
a decision problem, is to measure, say, the first qubit, 
and decide the problem according to this measurement 
outcome. In Ref. |55j it was shown that for every such 
decision problem, there exists another quantum circuit 
U' such that the evaluation of (0|?7'|0) is equivalent to 
the decision problem solved by applying U and measuring 
the first qubit. In this sense we have: 

Lemma 5 Knowledge of {0\U{il)\0) suffices to deter- 
mine the output of a quantum circuit which is being used 
to solve a decision problem. 

For a proof see, e.g., Ref. [55] . 



D. Proof of Theorem [T] 

Collecting everything we now prove our main theorem. 
We first need one more technical Lemma: 

Lemma 6 A quadratic form x'^Ax over GF(2) is linear 
in X (equal to a;*diag(^)) iff A is symmetric. 

Here diag(A) denotes a vector comprising the diagonal 
of A. The proof of this Lemma is presented in Appendix 

m 

Proof of Theorem [ij We start from the second claim 
of the Theorem, namely we prove that we can efficiently 
decide whether a given graph belongs to the set Q, and 
that we can find some w if it does belong. Let G be 
a given graph and let K be the matrix whose columns 
are a basis of Ker(yl), where A is the incidence matrix 
of G. This means that any a € Ker(^) may be written 
as a = Kx where x is an arbitrary m = dim(Ker(yl))- 
dimensional binary vector. Using this wc may rewrite 
Eq. (l34f over GF(2) as 



x*K^QKx = {Kxfi 



(36) 



Since the right-hand side is linear in x for all x, the 
left-hand side must also be linear in x. It follows 
by Lemma [6] that K*QK is symmetric, and more- 
over that the quadratic form x* K^QKx can be writ- 
ten as x*dia,g{K''QK). Thus, solving Eq. (34 1 for 
w is equivalent to solving the linear system x^K^w + 
x*diiig{K*QK) = 0, or 



x\K^w + diag(X*Qif)) ^ 0. 



(37) 



Since this equation must be true for all x, it follows that 
K'^w + dia,g{K^QK) = and hence that w is the solution 
to 



Since K is efficiently construct able, w can also be found 
efficiently using standard methods for solving linear equa- 
tions over GF(2). 

Now for the first claim of the theorem, which states 
that the circuits corresponding to the graphs in Q are 
CES. O is a downwardly closed set of graphs which gen- 
erates a set of compatible quantum circuits via the 
mapping described above. In turn we have the corre- 
sponding Ising model instances /(fieu;) and by assump- 
tion, the Gfl- compatible bond distributions w for the cir- 
cuits e rieiii and graphs G E Q. Lemma [l] states that 
these instances are planar. It follows that they are CES, 
i.e., we may compute the Ising partition function for any 
of these instances efficiently with a classical computer. 
This is due to a result by Kasteleyn, who gave a classical 
algorithm for the exact evaluation of the Ising partition 
function of any planar graph in the absence of an exter- 
nal magnetic field [58|. According to our definition of 
CES quantum circuits (Definition 12 1, all we need is to 



be able to obtain the evaluation in a time polynomial in 
the number of qubits (which translates to the number of 
vertices), and the desired number of bits of precision of 
Za(A), which is achieved by the algorithm given in |58j . 
Now, since {0\U{n)\0) cx Za(A) [Eq. and for any 



graph in Q we have an efficient way of classically deter- 
mining Za(A), we are thus able to determine the matrix 
element (0|[/(ri)|0) for any quantum circuit in fie^ ef- 
ficiently. It follows from Lemma [5] that any quantum 
circuit in ^Iqw which solves a decision problem is CES. 
■ 

We remark that this technique can be used to prove 
that quantum circuits which correspond to non-planar 
classes of graphs for which the Ising partition function 
has efficient classical evaluation schemes, e.g., graphs of 
bounded tree width, are CES. We suspect that some of 
the results obtained in Ref. [39' may be reproduced in 
this way. 

We further remark that due to algorithms for planarity 
testing, given a quantum circuit one can test if it belongs 
to the class of CES quantum circuits. For example, 
a simple test follows from the Eulerian criterion of pla- 
narity: \E\ < 3\V\ — 6 where \E\ is the number of edges 
and \V\ is the number of vertices. (This follows from 
the application of a handshaking lemma to the famous 
relation |i^| — -I- = 2, where F is the number of 
faces [57] . ) Examining the close relationship between the 
circuit representation H and the incidence matrix GH, 
one can give the restriction 

number of gates < 3 (number of qubits) — 6, (39) 

provided that the universal gate set consists of rotations 
about products of Pauli operations. A circuit f2 for which 
Eq. ( 39 1 holds generates a planar graph G via the map- 
ping we have described. Provided Eq. (34 1 has a non- 



K*w = diiigiK^QK). 



(38) 



trivial solution w, it follows that ft is Gw-compatible, 
and hence CES by planarity of G. 
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VI. FURTHER CHARACTERIZATION OF THE 
CLASS OF CES QUANTUM CIRCUITS Qe^ 

Our motivation in this subsection is to present a re- 
sult on CES quantum circuits which allows a compari- 
son to the recent results presented in [IS] and [3^. In 
both papers, results dependent on quantum gates being 
restricted to nearest neighbor qubit operations in one di- 
mension are presented. Via our construction we derive a 
similar result, but show that the restriction to nearest- 
neighbor operations and one dimension can be lifted. We 
begin with a simple example. 



A. Graphs in O are compatible with CES circuits 
which include non-nearest-neighbor operations 

Recent work in Ref. demonstrates that any cir- 
cuit that is built out of X-rotations and nearest neighbor 
Z(E)Z rotations can be efficiently simulated. Now assume 
that one is restricted to a class of planar graphs, 6p, for 
which the number of even subgraphs scales polynomially 
with the number of vertices. This restriction is not nec- 
essary and is introduced merely for simplicity. Let us 
call the corresponding set of quantum circuits (under the 
mapping presented above) flp. Upon inspection of the 
incidence matrix of a typical graph in Op one sees that 
even though the majority of incident vertices are nearest 
neighbor, there are several that are not, no matter how 
one labels the vertices. For example, consider the graph 
depicted in Fig. [T] 




FIG. 1: A graph with only one cycle. 



The incidence matrix is given by 

" 1 " 

1 1 1 

110 

1110 

1 

. 1 1 . 

and one possible circuit representation H is given by 

" 1 ' 

1 

1 

1 1 1 
110 
110 
1 
1110 

1 
1 
1 1_ 

Note that the fifth and sixth columns correspond to 
gates of the form 

respectively. The superscripts indicate which qubit is be- 
ing operated on, and thus one can clearly see that non- 
nearest neighbor interactions are possible. This example 
demonstrates that our construction may extend the re- 
sults in [42] , for example, by linking together graphs like 
the "necklace" shown. Note, however, that the nearest- 
neighbor restriction can only be relieved slightly as most 
of the interactions will in fact remain nearest-neighbor. 
Taking this example as motivation, what follows is a more 
general construction which will be used to pursue a better 
understanding of CES circuits. 

B. A class of CES circuits with non-nearest 
neighbor gates 

We now define a simple subclass of planar graphs which 
have a polynomial (in the number of vertices) number of 
even subgraphs. Because of planarity this class of graphs 
correspond to CES quantum circuits. However, note that 
this class of graphs is by no means an exhaustive charac- 
terization of all planar graphs which have a polynomial 
number of even subgraphs. 

Definition 13 A basis of the null space of the incidence 
matrix CH is referred to as a cycle basis. 

Definition 14 Let Opc be those planar graphs with V 
vertices and E = V + 0(log V'') edges, where k G M+. 

Proposition 1 Qp^ has a polynomial, in V , number of 
even subgraphs. 
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Proof. A cycle basis consists of a set of connected even 
subgraphs of a given graph in Qpc- The dimension of the 
null space (or the number of elements of the cycle basis) 
is in this case equal to the number of edges minus the 
rank of CH. Recall also that the rank of the incidence 
matrix equals the number of edges minus the number of 
components (which is one in our case). Thus, asymptot- 
ically, one has, 

nullity = V" + 0(fclogF)-rank(Ci/) = OilogV''). (40) 

Now, note that the null space allows all possible sums of 
the basis and therefore we are left with 0{V'') elements 
as claimed. ■ 

One could imagine graphs in Qpc as being sparse 
graphs consisting of cycles (even subgraphs) strung to- 
gether along trees without too many branching points. 
This is due to the relationship between E and V given 
above. One can see that a branch without a cycle always 
adds an additional vertex (one edge has two vertices) and 
the only way that the relationship between vertices and 
edges can be satisfied is if the number of branches is kept 
smaller than the number of cycles. That is, there will 
need to be more cycles than edges that do not terminate 
at a cycle. Further, the incidence matrix of these struc- 
tures, like the example above in Fig.[T] will have columns 
that consist of nearest neighbor consecutive "I's" for the 
majority of positions. This rule is broken when a tree 
branches and when one runs into a cycle. By Theorem 
[l] this can only happen 0{logV'') times. As E is the 
number of gates and V is the number of qubits in the 
corresponding quantum circuit, we have just proven: 

Corollary 1 A quantum circuit consisting of gates of the 
form e~*^('^* ^(s>y'-^>) ^ which act on nearest neighbor qubits 
except /or 0(fc log(#of qubits)) gates, which can act on 
qubits i and j such that |i — j | > 2, is CES. 

Note: Z operations may be included in the exponent 
of this operator at any position not occupied by the X 
and Y operations. 

The important point in the last corollary is that we 
allow non-nearest neighbor gates, thus extending the re- 
sults of [35], and also [M], where only nearest-neighbor 
CES quantum circuits were considered. 



C. O includes some but not all outer planar graphs 

So far we have stressed the special role of A'3 3 and i^s, 
which led us to the conclusion that all graphs in O are 
planar. However, it turns out that we can be more spe- 
cific, since we have also been able to show that Ki^ and 
K^_3 are forbidden minors for O (Lemmapl see Appendix 
[Cjfor a description of the proof, using mathematical soft- 
ware ). I n other words, there does not exist a solution to 
Eq. (34) for the graphs and ^"3^3. These graphs play 



Definition 15 For any planar graph, there are regions 
bounded by the cycles of the graph and an unbounded re- 
gion outside of all the cycles. An outerplanar graph is 
a planar graph for which every vertex is within the un- 
bounded region when it is embedded in the plane such that 
no edqes intersect. 



a role in characterizing the set of outerplanar graphs |57j , 
which we define next: 



For example, the graph in Fig. [T] is outerplanar. More 
informally, a graph is outerplanar if it can be embedded 
in the plane such that all vertices lie on the outer (exte- 
rior) face. A graph G is outerplanar iff Ki + G (a new 
vertex is connected to all vertices of G) is planar [BIT . The 
characterization of relevance to us is the following analog 
of Kuratowski's theorem for planar graphs (described in 
Appendix |A| : 

Theorem 2 (Chartrand &; Harrary [61]) A graph 
is outerplanar if and only if it has no subgraph homeo- 
morphic to or K2.3. 

In other words, if 2, 3 is a forbidden minor for outerpla- 
nar graphs, where K2^3 is like 3 except that one side 
of the bipartite graph has two vertices instead of three. 

Proposition 2 If a graph is outerplanar then it does not 
have -H'3,3 as a minor. 

Proof. Assume that ^3^3 is a minor of some G G 
Outerplanar. Then G has ^^2,3 as a minor, since -ftr2.3 
is a minor of ^3,3. (This is easy to see: just contract 
one edge and delete another.) But by Theorem [2] outer- 
planar graphs cannot have ii'2,3 as a minor, which is a 
contradiction. Thus G cannot be outerplanar. ■ 

Note that the converse is not necessarily true, i.e., not 
all graphs which do not have -^^3.3 as a minor are outer- 
planar. 

Proposition 3 Outerplanar Graphs n O 7^ and 
Outerplanar Graphs 7^ O. 

Proof. By Lemma|2] if G G O then G cannot have K4 or 
3 as minors. By Theorem [2] and Proposition |2] if G' 
is an outerplanar graph then it cannot have K4 or -^'3^3 
as minors either. This suggests that O may have graphs 
in common with the set of outerplanar graphs. We have 
verified, using mathematical software, that the intersec- 
tion is indeed nonempty. For example, we have found 
that certain trees with cycles, which are outerplanar by 
construction, are in @. Moreover, we have verified using 
mathematical software that ^2,3 € O, i.e., it is not a for- 
bidden minor for the existence of a solution w. But -fir2^3 
is not outerplanar, hence there are graphs in Q which are 
not outerplanar (in particular, all subdivisions of K2.3). 
m 

This is interesting since some problems are NP- 
complete for subclasses of planar graphs but solvable in 
polynomial time for outerplanar graphs. Some exam- 
ples are the chromatic number, Hamiltonian path and 
Hamiltonian circuit. Another example is the page num- 
ber which is one for outerplanar graphs. This means 
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that we can embed the vertices on a Hne which divides 
the plane into two subplanes, and draw all edges in one of 
the subplanes without crossing. In a sense, outerplanar 
graphs are "easy" computationally. The fact that Q in- 
cludes non-outerplanar graphs thus suggests that it may 
include interesting computational problems. 



seems likely that also non-planar graphs (i.e., those con- 
taining or ^^3,3 as minors) may have a satisfying w 
according to Eq. (42 1. We leave this as a problem for 



future investigation. An important point concerning the 
more general condition ( 42 1 is that wc do not know if it 



can be efficiently tested for a satisfying bond distribution 
w, let alone solved for such a w. 



VII. DISCUSSION AND FUTURE DIRECTIONS 



In this section wc briefly discuss two possible future 
directions for research. 



A. General condition for the bond distribution 



So far we have assumed the sufficient condition (|34| in 
order to obtain the desired equality between the parti- 
tion function and the circuit matrix element (0|J7(ri)|0) 
[Eqs. (25 1 and (32 1]. We have also shown that a satisfy- 



ing bond distribution w can always be efficiently tested 
for an computed, under Eq. (34 1. Let us now relax the 



constraint of Eq. ( 34 1 by considering a more general way 



in which the desired proportionality. 



aGkcr A 



aekcr CH 



(41) 

can be obtained. Indeed, Eq. (p54| is clearly not a neces- 



sary condition. The following construction demonstrates 
that it is likely that the number of cases which do not 
have a solution w for the bond distribution is much 



smaller than the case we analyzed given by Eq. ( 34 1 



Note that in Eq. (41) the powers of the A's are the 



weights of the null vectors a, that is the number of ones 
in a. Thus it is possible for an equality to occur for a 
given term in the sum in the two sides of Eq. ( 41 1 for 



different a's, as long as the weights of the a's are equal. 
This gives us the constraint for the following. One can 
organize all the a's in bins in terms of weights from 1 to 
\E\ = N. Let us now take bin r, i.e., the set of vectors of 
weight r. Let a^i , . • . , ctm be all the null vectors of CH 
of weight r. Then, if for all r 



{a*j^lwtr(iJ'CiJ)arji = a,.j2 • w} A 
{a^j-^lwtr(iJ*CiJ)arj2 = a»-j3 ■ w} A ■ 
{a* ^lwtr(i/*CiJ)arj„ = arj^ ■ w}, 



(42) 



where {ii,-.-,jn} is any permutation of the numbers 
{l,...,n}, then Eq. (41 1 would be satisfied. Clearly, 
Eq. (34 1 is a special case of this more general condition. 



This demonstrates that it is likely that a satisfying bond 
distribution w can be found for a given graph even if 
the sufficient condition (34) cannot be satisfied. Loosely, 



this is due the fact that there are many conjunctive state- 
ments of the form ( 42 ) that can be satisfying. In fact it 



B. Computing the Ising partition function 



As mentioned in Section IV A Eq. (35) has two con- 



sequences, and our focus in this paper has been on the 
ability to find CES circuits using known results about the 
hardness of computing partition functions. Let us now 
briefly consider the other consequence, namely the fact 
that if we are able to determine (0|[/(il)|0), then we are 
able to determine the partition function Z/^(X). 

A fully-polynomial randomized approximation scheme 
(fpras) for the fully-ferromagnetic Ising partition func- 
tion was presented in |62j . It is well known that having 
an fpras for the non-ferromagnetic Ising model implies 
that NP = RP (randomized polynomial time) which 
would be quite unexpected [SU] . It should therefore be of 
no surprise that no fpras for this problem has been found, 
even with quantum resources. However additive approx- 
imation schemes seem likely and in fact one was given in 
pS] for the related Potts model partition function, even 
though the instances that they were able to account for 
are not known to be BQP-complete and the hardness is 
in fact unknown. 



Equation ( 35 1 precisely relates a matrix element of a 



quantum circuit with the value of the partition func- 
tion of the Ising model for a corresponding graph in- 
stance. This means that if we could approximate the 
matrix element, we would have an approximation for the 
Ising partition function. Due to the Hadamard test, it 
is well known that a polynomial estimation of this ma- 
trix element is BQP-complete. (See Ref. |53j for a de- 
scription of the Hadamard test.) Specifically, by making 
1/e^ measurements, one can either have Re(0|J7(O)|0) or 
Im(0|C/(ri)|0) to precision e, but one must keep in mind 
that this approximation is an additive one. This means 
that with some probability of success bounded below (say 
by .75) the approximation returns m such that 

(o\uin)\o)-s-p<m< {o\u{n)\o) + s -p, (43) 

where p is a polynomially small parameter and S is 
the approximation scale of the problem. Note that if 
6 — O((0|J7(ri)|0)) then the approximation will be an 
fpras ^23j . Equations ( 35 1 and (E3|) taken together quan- 



tify how a measurement of (0|t7(T2)|0) yields an approx- 
imation of the partition function of the Ising model in- 
stance A corresponding to the circuit fi. 
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VIII. CONCLUSIONS 



APPENDIX A: ESSENTIAL ELEMENTS FROM 
GRAPH THEORY 



We have provided a construction that allows one to 
determine if a given quantum circuit corresponds to a 
class of quantum circuits which are classically efficiently 
simulatable (CES). This was done by looking at the cor- 
responding graph instances of the classical Ising model 
using a mapping previously introduced in [2? . This was 
then used to conclude that any class of quantum circuits 
which solve decision problems and are restricted to cer- 
tain planar graph instances are CES. Our main result is 
stated in Theorem[l] which characterizes the class of CES 
circuits via the set of planar graphs O. We have given a 
partial characterization of by stating that its obstruc- 
tion set includes ^3,3 and K4 (and hence, by downward 
closure also ^^3,3 and K^). An interesting open problem 
is to give a complete characterization of the obstruction 
set; we know from the Robertson Seymour Theorem that 
this set is finite, since we have proved that Q is down- 
wardly closed. 

Our mapping can also be used to construct a quantum 
algorithm for the additive approximation of the partition 
function. However, there are two issues. The instances 
we are able to handle are constrained by our use of equa- 



tion ( 34 1 which does not capture all the ways a certain 



bond distribution may satisfy equation (41), but which 



simplifies our analysis greatly. The other issue is the fact 
that our mapping may fail to provide information about 
the bond distribution of a given graph. 

An open problem is to obtain a better understanding 
of what the complexity of finding the bond distribution 
for a particular graph instance is. This understanding 
will have consequences in our knowledge of where BQP 
is in the complexity hierarchy, as we will be able to relate 
the simulat ability of universal quantum circuit families 
with the complexity of finding bond distributions. On the 
other hand, it is possible that the complexity of finding 
bond distributions is somehow incorporated in the power 
of the quantum circuit that corresponds to the graph 
instance of the Ising model, in the sense that the circuit 
corresponding to a planar graph (under our mapping) 
may not be CES, because the effort of obtaining the bond 
distribution via Eq. (42) blocks such a simulation. 



Here we review some essential definitions and theorems 
form graph theory needed for the results presented in 
this work. A good reference for these concepts is the 
wikipedia article on planar graphs, or Ref. |57j . 

Definition 16 A subgraph H of a given graph G is called 
a minor (or child) of G if it is isomorphic to a graph 
that can be obtained from G via a sequence of edge dele- 
tions, edge contractions, or deletion of isolated vertices. 
Edge contraction is the process of removing an edge and 
combining its two endpoints into a single vertex. Edge 
deletion removes an edge without removing its vertices. 

Definition 17 The set of graphs S is downwardly closed 
with respect to minor ordering if whenever G is a member 
of S , then so is any minor of G. 

A trivial consequence of the definition of downwardly 
closed sets is that every such set has an obstruction set: 

Corollary 2 An obstruction set for G is a set of minors 
of G, also called forbidden minors, with the property that 
they prevent downward closure. 

In other words, if one constructs a set of minors of 
G and encounters a minor that violates the property for 
which downward closure is being tested, such a minor is 
called forbidden, and belongs to the obstruction set. Now 
comes a seminal theorem due to Robertson and Seymour 

m- 

Theorem 3 (Robertson-Seymour) Every downwardly 
closed set of graphs (possibly infinite) has a finite ob- 
struction set, i.e., a finite set of forbidden minors. 

For our purposes an important example are the planar 
graphs: 

Definition 18 A planar graph is a graph which can be 
embedded in the plane, i.e., it can be drawn on the plane 
in such a way that its edges may intersect only at their 
endpoints, i.e., edges never cross. A nonplanar graph is 
a graph which cannot be drawn in the plane without edge 
intersections. 



Acknowledgments 



This material is based upon work supported by the Na- 
tional Science Foundation (NSF) under Grant No. PHY- 
0802678, and by the Army Research Office under grant 
W911NF-05- 1-0440. DAL thanks the Institute for Quan- 
tum Information (IQI) at Caltcch where part of this work 
was done. IQI is supported by the NSF under Grant No. 
PHY-0803371. 



There are two particularly important nonplanar 
graphs, denoted (the complete graph on five vertices; 
complete means that each pair of vertices are connected 
by an edge) and K^^ (the complete bipartite graph on 
six vertices, three of which connect to each of the other 
three). They are depicted in Fig. [2] 

Planarity is characterized by Wagner's Theorem 57J: 

Theorem 4 (Wagner) A finite graph is planar if and 
only if it does not have or 3 as a minor. 
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FIG. 2: The various graphs playing a role in defining the 
obstruction set for 6 and outerplanar graphs. 



In other words, if one deletes or contracts the edges of 
a graph and finds one of these minors, then the graph 
is not planar. Hence -fCs and ^ are forbidden minors 
(form an obstruction set) for planarity. 

For completeness we note that an alternative charac- 
terization can be given in terms of the concept of a sub- 
division of a graph: 

Definition 19 A subdivision of a graph results from in- 
serting vertices into 



Thus, while deletion and contraction of edges shrinks 
a graph down, subdivision builds it up. 

Definition 20 Two graphs G and G' are homcomorphic 
if there is an isomorphism from some subdivision of G to 
some subdivision of G' . 

Theorem 5 (Kuratowski) A finite graph is planar if 
and only if it has no subgraph homeomorphic to Kr, or 
K3.3. 



In other words, a finite graph is planar if and only if 
it does not contain a subgraph that is isomorphic to a 
subdivision of or ^^3^3. 

In this work we give a criterion for graph member- 
ship (in a certain set Q which is defined in Definition 11 ) 



based on the existence of a solution for the system of lin- 
ear equations over GF{2) defined by Eq. (34 1. We want 
to be guaranteed that this membership has an ordering 
in the sense that if G is a member, then so is every minor 
of G. In other words, we are testing for downward clo- 
sure. The Robertson-Seymour theorem guarantees that 
membership in our set Q is not obstructed by an infinite 
set of graphs. However, the situation is in fact far better: 
membership of a graph to a fixed downward closed set can 
be checked by running a polynomial time algorithm for 
all elements of the obstruction set (if it is known), since 
searching for a minor on a given graph only requires cu- 
bic time |57j . In fact, checking whether a graph is planar 
can be done in linear time. 



For completeness we include the definition of a hyper- 
graph as our correspondence between quantum circuits 
and graphs is actually a mapping between circuits and 
hypergraphs pS] . 



Definition 21 A hypergraph is a generalization of a 
graph where edges are replaced by hyperedges. Let 
V — {vi, V2, ■ ■ ■ ,Vk} be the set of vertices and let E = 
{ei, 62, . . . , e„} be the set of hyperedges. Each Ci — 
{^iiiVi2, ■ ■ ■ ,Vim} is a collection of vertices where each 
v^., G V. 



Thus the main difference from ordinary graphs is that 
edges consist of arbitrary collections of vertices rather 
than two and thus graphs are special cases of hyper- 
graphs. As shown in Ref. [IH], the existence of hyper- 
edges is what gives us access to the universal gate set 
presented in [49^ , via the two assumptions enumerated at 
the end of Section |IV[ However, the circuits that corre- 
spond to graphs is what is interesting here, as our results 
depend on information about the Ising partition function 
defined on ordinary graphs. 



APPENDIX B: PROOF OF LEMMAS 

Here we present the proof of lemmas [T] [3] and [6] which 
we repeat for convenience: 

Lemma [S] // a graph G is a member of Q, then so is 
G \ Cj or G/cj, i.e., the deletion or contraction of an 
arbitrary edge e.j from a graph in O is also in Q. 
Proof. Assume that G G O. Recall that a graph G is an 
element of O if there exists some solution w to the set of 
linear equations over GE{2) 



w = a 



(G) 



(Bl) 



where A^*^' is the matrix whose rows are elements, a^, 
of the nullspace of the incidence matrix of the graph G 
(given as the Ising instance), and a^*^-* is the vector whose 
entries are the a*lwtr(i/*Gi/)ai. These null elements, Oi, 
correspond to the even subgraphs of G and will be re- 
ferred to as cycles (recall Definition [4]) . From elementary 
linear algebra we know that a solution exists if a'-'^^ may 
be written as a linear combination of columns of A^'^^ or 
in other words if 



Rank[y4(^)|a^^^] = Rank[A'-^^]. 



(B2) 



We must demonstrate that after we either delete or con- 
tract an edge, and arrive at the subgraph G', we have 



Rank[^(^'V^'^'^] - Ra.nk[A^^">]. 



(B3) 
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We shall demonstrate this with an edge deletion as the case of a contraction is similar. We begin with a given graph 
G which is a member of F^. We have 



(G"*CG')^'^^ = 



Gf r^l 

11 '-'21 

^12 ^22 

Is ^2s 



Gv2 

c 

^ vs 

r" 



^21 


^22 ■ 


■ ■ 


r" 





• 


• • 


■ 




^42 ■ 


• G'^g ■ 


• G'^jy 





• 


■ • 


■ 







r" 





g: 



vN 





(B4) 



Using Einstein notation this equals 



2i-l,l"-^ "-^2i-l,l"-^ 



,2i,2 



2i-1.2^ 



^2i-\,N^ 



2i-\,2'^ 



^2i-\.N^ 



2i-l,l"-^ 

Gf /~it'- 
2i~\,2^ 



'^2i-l.N'^ 



2i,s 



Gl /^l 
2i-l,l"-^ 

^'21-1,2^'' 



2i,N 
2i,N 



'-^2i-l.iV^' 



2i.N 



(B5) 



Keep in mind that if we were to delete an edge from G, t his w ould correspond to losing a column from CH which 
would correspond to losing, say the s"^ column from matrix (B5l. Taking the lower triangular portion of matrix (B5l 
and calculating we find that the m**^ element of the vector a^'~'> is 



0-m2[o-mlG'2i_i2 



G' 



2i,ln 



a. 



m3 



[a„iiG2j_i 3G' 



2j,l 



/2i,l 



a. 



m{k 



-1)^21-1 hG' 



0'm2G'2i-i^zG' 

y2i,fc-li 



2i,2-| 



-i,2i,JV-li 



where ami is the i"^ element of the m 
Now, let 



amn[amlG2j_ljV^' ' + ' ' ' + '^m(Ar-l)G"2j_i jy^' 

null vector oi GH or the matrix element 



(B6) 



a. 



m(s+l)'2msG2j_x^s4. 



-iG' 



2i,s 



^m(s+2)ttmsG'2i_i ,5+2 



,G' 



2i,s 



' ^mN^msG2 



2i-l,N 



G' 



2i,s 



(B7) 



(G) 

This is the portion of am that would vanish if we were 
to omit the edge that corresponds to the s*"^ column of 
the matrix (B5 1. This means that if we remove this edge, 
we will end up with the subgraph G' and we can write 

«l?^-«if)+a. (B8) 

This equation is saying that the m**^ entry of the right 
hand side of 



w 



(B9) 



is given by the m entry of the right hand side of 

= (BIO) 

(the corresponding system of equations for the graph G') 
plus the term 

From the assumption that G G F^ and by construction 
we have 



JG') 



(G') 



£1 



where the 6i are coeflicients in GF{2) and the Ci are 
columns of A^^), i.e., Rank[A("^) Rank[A('^)]. 
Thus we have 



a 



(C) 



(B12) 



How does the matrix A change as we go from G — > G' 
by this edge deletion? If the edge is a dangling edge, i.e., 
not part of a cycle, then we lose a column (column s) but 
if the edge deletion causes the breaking of M cycles, then 
A will lose M rows (in addition to column s), as the rows 
encode the cycle structure of the graph. In this case, the 
dimension (or length) of a^'^ ^ will be M less than the 
dimension of a^'^^ and the Ci will also be shorter by M 
entries. We call these shorter Ci, c[. Further, and most 
importantly, £* will vanish, as mentioned. After taking 
this into consideration we now can conclude that 



a 



(C) 



(B13) 



(BU) where = [( 



1^2 



■c'jy_i]- Thus, 



Rank[^(^'V^'^'^] = Ra.nk[A^^">]. (B14) 
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The proof for edge contractions is similar. The main 
difference is that an edge contraction does not cause the 
loss of a cycle except when the edge in question belongs to 
a cycle of length three. Thus in general, the contraction 
case is simpler except when dealing with cycles of length 
three. In this case the proof caries over in the same way. 
■ 

Thus the set of graphs Q is downwardly closed. 
Lemma [l] The obstruction set for Q is finite. 
Proof. The set of graphs Q is downwardly closed by 
the above lemma. One may then apply the Robertson- 
Seymour Theorem (Theorem [s]) and immediately con- 
clude that the number of forbidden minors of O is finite. 
■ 

Lemma [b] A quadratic form Ax over GF(2) is linear 
in X (equal to a;*diag(A)) iff A is symmetric. 
Proof. Let x be an m-dimensional column vector and 
A an m X m matrix, both over GF(2). Consider the 
quadratic form x^Ax — ^j^- XiAijXj and assume that A 
is symmetric: A = A* . Then 



i i>j 

, AijXiXj -)- ^ ^ AiiXi ~\- ^ ^ AijXjxiJiXh) 



E 



j>i 



where in the second line we used xf — Xi [true over 
GF(2)], exchanged i and j in the third summand and 
used Aij = Aji. The first and third summands are equal 
and hence add up to zero over GF(2). We are left with 
the, second, linear term, i.e., x^Ax = a;*diag(A), where 
diag(A) denotes a vector comprising the diagonal of A. 

Next, assume that A is not symmetric. Then there 
exists a pair of indices i' < j' such that Aiiji ^ Aj'ii, 
i.e., Aiiji + Ajui = 1. As above, we have: x^Ax = 
Y.i<j AijXiXj + ^ii^i + Hj>i AjiXjXi. Consider the 
index pair (i', j') in this sum: 

Ai'j'X^'Xj' Aj' i' Xj' Xi' (^A^'j' ~{- Aj'i'^Xi'Xj' x^'Xj'. 

(B16) 

Thus the quadratic form contains at least one non-linear 
(quadratic) term Xi'Xji. ■ 



APPENDIX C: ALGORITHM FOR MINOR 
TESTING 

Note that all calculations are done modulo 2. 

Input: A graph G for which we wish to determine 
if there exists some satisfying edge interaction w that 
satisfies Eq. (34 1. Specifically, we considered G — i^3,3, 
K5, 3 (-K^3,3 with one edge deleted) and K4. 



Output: A binary vector w that is a satisfying bond 
distribution or a null vector (indicating no such bond 
distribution) . 



1. From the incidence matrix A of G obtain the fol- 
lowing items: 

(a) All vectors belonging to the null space £ of 
the incidence matrix. These row vectors form 
a matrix M. 

(b) Construct a matrix representation H of the 
possible corresponding quantum circuits (un- 
der the mapping presented earlier). From H 
construct Q = IwtiH^CH. This matrix will 
have variables Zi corresponding to all the pos- 
sible ways that one can include or omit Z op- 
erations (changing these affects the types of 
edge interactions that one obtains, if any.) 

2. Form the vector B whose i*'' entry is Bi — a\Qai, 
where are the elements of the null space of C. 
This is the left-hand side of Eq. ( 34 1 . Each entry 



of B consists of linear equations whose variables Zi 
represent the presence or absence of a Z operation 
in a quantum circuit that corresponds to G. 

3. Form a matrix W whose rows are all possible bonds. 

4. Produce a matrix D whose i'^ row Di is equal to 
MWi. These are all possible values of the right 
hand side of Eq. ( p4| . (Note that due to symmetry 
there will be many repeats, so that the total num- 
ber of possible bonds to check is far fewer than all 
possible bonds.) 

5. Attempt to solve the system of linear equations 
Bk =Di^k,Bk= D2,k, ...,Bu = D\Lik over GF(2) 
(where k runs from 1 to the number of rows of D) 
for the variables Zi. A solution for some k gives 
information for a specific circuit representation H . 
If no solution for the Zi exists, then there is no sat- 
isfying bond distribution w for Eq. ( [34| . This is 
indeed the case for if 3^3, K^, Kj, ^, and K^. If there 
is a solution for some fixed k, continue. 

6. Take this specific H (i.e., this H has no vari- 
ables and corresponds to a specific circuit given 
by the solution of the Zi above) and now form 
Bi = al{lwti H*GH)ai again. This time however, 
B contains no variables and is a numerical vector. 
Thus, one now has the binary vector B and the 
matrix M. 



7. Solve the linear system AIw 
edge interaction w. 



B and output the 



Having applied this algorithm we proved using math- 
ematical software that the obstruction set for O includes 
7^3 3, K5, ^3,3 and K4, i.e., Lemma[2j 
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